Flux qubits with neutral currents in optical lattices 
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We study an experimentally feasible qubit system employing neutral currents. Our system is based on bosonic 
cold atoms trapped in ring-shaped optical lattice potentials. The lattice makes the system strictly one dimen- 
sional and it provides the infrastructure to realize a tunable ring-ring interaction. By breaking the Galilean 
invariance we demonstrate how atomic currents through the lattice provide a realization of a qubit. We break 
the Galilean invariance either by artificially creating a phase slip in a single ring, or by considering two ho- 
mogeneous ring lattices, coupled by tunneling interaction. The Hamiltonian of the system effectively leads to 
a washboard potential in the phase representation, tilted by the applied 'flux'. The single qubit infrastructure 
is experimentally investigated with tailored optical potentials. An experimentally feasible scheme of the two- 
ring-qubit is discussed. In this case, the dynamics is demonstrated to show macroscopic quantum self trapping. 
Time-of-flight expansion maps the pattern of atomic currents into detectable atomic density distributions. Based 
on our analysis, we provide viable protocols to initialize, address, and read-out the qubit. 



A qubit is a two state quantum system that can be ma- 
nipulated coherently despite the environment, coupled to its 
neighbors, and measured. All proposed experimental im- 
plementations present qubit specific advantages and difficul- 
ties. These systems include ion traps[l], NMR[2], Joseph- 
son junctions [3], quantum dots [4] and cold atoms [5, 6]. In 
particular cold atoms and solid state Josephson junctions im- 
plementations have complementary characteristics. In neu- 
tral cold atoms proposals the qubit is encoded into well iso- 
lated internal atomic states. This allows long coherence times, 
precise state readout and in principle scalable quantum regis- 
ters. However individual qubit (atom) addressing is a deli- 
cate point for optical lattice architectures [7]. Qubits based 
on Josephson junctions implement the two level system either 
using charge or by clockwise-anticlockwise superconducting 
currents. Such devices allow fast gate operations and make 
use of the precision reached by lithography techniques for cir- 
cuit design [8]. The decoherence, however, is fast in these 
systems and it is experimentally challenging to reduce it. For 
charge qubits the main problem arises from dephasing due to 
background charges in the substrate[9]. Flux qubits are insen- 
sitive to the latter decoherence source, but are influenced by 
magnetic flux fluctuations due to impaired spins proximal to 
the device [10]. 

Here we aim at bringing together the advantages char- 
acterizing Josephson junctions and cold atoms. The ba- 
sic idea is to use the neutral currents flowing through ring 



shaped optical lattices to realize the cold atoms analog of 
flux qubits[ll]. Persistent currents in optical ring lattices are 
experimentally generated in different ways: by rotating the 
condensates [12, 13], by exploiting dark- states [14- 16], or by 
implementing the phase slips as a Berry phase[ll, 17]. Su- 
perpositions of persistent currents have been thoroughly in- 
vestigated in a series of papers by Rey et. al In the case of 
homogeneous rings the superposition state can in principle oc- 
cur but with important restrictions on the density of superfluid 
and strength of the interaction; further, the visibility of the two 
states vanishes in the thermodynamic limit[18, 19]. Some of 
these limitations can be overcome by breaking the Galilean 
invariance of the system[20, 21]. 

In this paper we demonstrate how persistent currents flow- 
ing in a ring shaped optical lattice [1 1] can provide a physical 
implementation of a qubit. The lattice potential plays an im- 
portant role in our approach because it provides the confine- 
ment of the atoms in a strictly one dimensional ring. Further 
it provides the means for precise control of the confinement 
and facilitates the qubit-qubit interaction. In our system we 
break the Galieian invariance. For a single ring this is realized 
by creating a localized 'defect' barrier along a homogeneous 
lattice [22]. Furthermore we prove that a qubit can be achieved 
with two homogeneous interacting rings arranged vertically 
on top of each other. In such a system the Galilean invariance 
is broken along the direction transverse to the two rings. For 
this scheme we analyze the real time dynamics and time-of- 



flight density distributions. 



Breaking the Galilean invariance on the single ring- We 
consider bosonic atoms loaded in a ring-shaped potential with 
identical wells, but with a dimple located at the site N - 1 
(see Fig.l), and pierced by a 'magnetic flux' $. The system is 
described by the Bose Hubbard Hamiltonian 
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BH 



^Y.ni{n,-l)- Y,Uie^'"''4^i.i + h.c.) (1) 



where a^'s are bosonic operators for atoms trapped in the the 
ring and n^ = alai. The parameters U describe the tunneling 
between the wells along the ring. Since the wells are all iden- 
tical but one, U = t^\fi = . . .N - 2 and tN-i = t' • Finally, 
U describes the s-wave scattering interaction [26]. The 'mag- 
netic flux' is <l> = /j^'^^ A(z)dz, where A(z) is the effective 
vector potential. The effect of the dimple is to induce a phase 
slip at the site N-l. We assume that the density of super- 
fluid is large enough to neglect the fluctuation of the number 
of atoms in each well. In this regime we can assume that the 
system dynamics is characterized by the phases of the super- 
fluid order parameter ^^'s, described by the quantum phase 
model [27] with Josephson coupling Ji = {n)ti ((n) is the av- 
erage number of bosons in each well). The magnetic flux $ 
can be gauged away everywhere but at the site (A^- l)-th [28- 
30]. Accordingly the phase difference along nearest neighbor 
sites can be considered small in the 'bulk' and the harmonic 
approximation can be applied. The partition function can be 
written as a path integral: Z = f V[(l)]e~^^^\ where the 6'[^] 
is the Euclidean action. Adapting from the approach pursued 
by Rastelli et al. [33], all the phases (j)i except 6 = (pN-i - 0o 
can be integrated out (the integrals are gaussian). The effec- 
tive action reads 
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2U{N-1) 

with the potential U(0) = (0 - <l>)^ - J' cos(i9). For 

large {N - 1)J'IJ and moderate N, U{9) defines a two- 
level system. The degeneracy point is <l> = tt: The 
two states are provided by the symmetric and antisymmet- 
ric combination of counter-circulating currents correspond- 
ing to the two minima of U{9). We observe that break- 
ing the Galilean invariance of the system provides an inde- 
pendent parameter J' facilitating the control of the poten- 
tial landscape. The interaction between and the (harmonic) 
bulk degrees of freedom provides the non local term with 
G{t) = Y.'^QY{uji)e'^^^^, uoi being Matsubara frequencies 
(iv-2)/2 l + cos[2^fe/(7V-l) 



and Y{uJi) 



k=l 



2J/7(l-cos[27r/c/(A/'-l)]) 



The external bath vanishes in the thermodynamic limit and the 
effective action reduces to the Caldeira-Leggett one [33, 34]. 
Finally it is worth noting that the case of a single junction 




FIG. 1. Realization of a ring-lattice potential with an adjustable 
weak link. We create the optical potential with a spatial light mod- 
ulator (SLM) which imprints a controlled phase onto a collimated 
laser beam. The figure above shows a measured intensity distribution 
(arb. units) with horizontal axes units in /im. The center peak is the 
residual zero-order diffraction. The structure is scalable and a lower 
limit is imposed by the optical wavelength. The ring lattice resem- 
bles the single qubit potential as discussed in Eq.(2). Calculation of 
the SLM phase pattern (kinoform) has been carried out using an im- 
proved version of the Mixed-Region- Amplitude-Freedom (MRAF) 
algorithm[22, 23] with angular spectrum propagator [24]. This allows 
us to calculate numerically the wavefront propagation without resort- 
ing to paraxial approximation. For calculation of the holograms an 
undistorted Gaussian wavefront as input intensity distribution for the 
kinoform has been assumed. A region outside the desired ring lattice 
pattern is dedicated to collect unwanted light contributions resulting 
from the MRAF algorithm's iterative optimization process. 



needs a specific approach but it can be demonstrated consis- 
tent with Eq.(2)[35] . 

Breaking the Galilean invariance with two homogenous 
coupled rings-We consider bosonic atoms loaded in two iden- 
tical rings Fig. 2. The system is described by Bose-Hubbard 
ladder: H = H^^ + Hj^^+Hint, where H^^ ^ are the Hamil- 
tonians as in Eq.(l) describing the bosons in the rings a and b 
respectively. The interaction is: 
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We observe that along each ring the phase slips imply twisted 
boundary conditions and therefore they can be localized to a 
specific site, say the N - 1-th. Following a similar procedure 
as employed above, the effective action reads 
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where each G^ (r ) is given by the expression found above for 
the case of a single ring. In this case the phase dynamics is 



provided by the potential 
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-Jcos[0a -0,- ^^(^a - ^6)] • (5) 

with J = {n)g. We observe that, for large N, the potential 
U{OajOb) provides that effective phase dynamics of Joseph- 
son junctions flux qubits realized by Mooji et al (large N's 
corresponds to large geometrical inductance of flux qubit de- 
vices) [32]. In there, the landscape was thoroughly analyzed. 
The qubit is made with superpositions of the two states \6i) 
and 1 6^2} corresponding to the minima of U{0a^0b). The de- 
generacy point is achieved by ^^ - ^a = tt. We comment 
that the ratio J/ J controls the relative size of the energy bar- 
riers between minima intra- and minima inter- 'unit cells' of 
the {Oa^O})) phase space, and therefore is important for de- 
signing the qubit. In our system J jj can be fine tuned with 
the scheme shown in Fig. 2. 

We now study the dynamics of two homogeneous rings cou- 
pled by tunneling. We will show that the density of the con- 
densate in the two rings can display characteristic oscillations 
in time. We make use of the mean field approximation to an- 
alyze the (real time) dynamics of the Bose-Hubbard ladder 
Eqs.(l), (3) (assuming that each ring is in a deep superfluid 
phase). Accordingly Gross-Pitaevskii equations are found for 
the variables (/:?a,i( -5 ) = {cii{s)) and(/P5,i(5) = (6^(5)}. Assum- 
ing that Oa - V^a,i+i -^a,i in each ring is site-independent, we 
obtain 

dz 

ds 

A + Xpz + ^ cos 9 (6) 



-V 1 - z^ sinO 



ds 



Vl^ 



where z = {Ni^ - N a)/ {N a -^ Ni^) is the normalized imbal- 
ance between the populations Na and Nij of the two rings, 
9 = Oa - 0}) and s = 2gs (h = 1). The parameters are 
A = [fia - l^b + t (cos ^- cos j^)]/g, A = U/{2g), and 
p = (Na + Ni))/N is the total bosonic density (we included 
the chemical potential /i^). Eqs.(6) can be solved analytically 
in terms of elliptic functions [36]. Accordingly, the dynamics 
displays distinct regimes (oscillating or exponential) as func- 
tion of the elliptic modulus k, depending in turn on A, A, and 
on the initial population imbalance 2:(0) = zq. Here we con- 
sider the dynamics at Xp « A, i.e. small U/g (the analy- 
sis of the solutions of the Eqs.(6) in different regimes will be 
presented elsewhere). The results are summarized in Fig. 3. 
We comment that, comparing with A = 0, the oscillations do 
not average to zero (therefore yielding a macroscopic quan- 
tum self trapping phenomenon [36]) and they are faster. 

The pattern of the circulating currents along the two cou- 
pled rings can be read out through the analysis of the time-of- 
flight density. As customarily, the spatial density distribution 
in the far field corresponds to the distribution in the momen- 
tum space at the time when the confinement potential is turned 
off: 
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FIG. 2. The proposed setup for the ring-ring coupling. Two paral- 
lel gaussian laser beams (G1,G2) with wavelength Ai are produced 
by a combination of two polarizing beamsplitter (BSl, BS2). The 
beam separation D can be controlled by moving mirror Mi . The 
beams pass through a lens with focal length / and interfere to form 
a lattice in z-direction. The distance between the lattice planes is a 
function of IjD [31]. The resulting one dimensional lattice is com- 
bined with two counter propagating Laguerre-Gaussian laser beams 
(LGl, LG2) of amplitude Eq and wavelength A2. To avoid inter- 
ference between the lattice in z-direction and the Laguerre-Gaussian 
beams, the wavelength difference between Ai and A2 should be large, 
typically a few nm. The inset shows the ring lattice potentials sepa- 
rated by d = Xif/D: Vutt = 4:Eo{fpi cos {kLoz)^ + cos (kcz)'^ + 
2 f pi cos {k lgz) cos (kcz) cos ((1)1)), where fpi are related to La- 
guerre functions [11]. In the figure / = 6 and p = 0. The ring-ring 
distance is adjustable by varying the distance D. The WKB estimate 
of the tunneling rate gives g ^ exp {-AE^fiXf jD). Such a system 
provides an effective two level system that can be exploited as a qubit 
(see text). 



P(k) = '"^^"^^'f,^'"^^^' E E E (7) 

cos[k|| -xii +{q+^){(t)i-(t)j)]{alaq) + 

cos [k|| • X|| + ((7 + ^)(0i - ^j)Wq^q) + 

2 cos [k|| . X|| + k,D + {q+ -^)(l)i -{q+ ^)0^)](46^) 



where w{kx^ky^kz) are Wannier functions (that we consid- 
ered identical for the two rings), ky • xy - kx{xi-Xj) + ky{yi- 
yj), Xi = cos 0i, ^i = sin (/)^ fix the positions of the ring wells 
in the three dimensional space, ^^ = 27ri/N being lattice sites 
along the rings. The density Eq.(8) is displayed in Fig.4 

Conclusions- We discussed a construction of flux qubits 
with atomic neutral currents flowing in ring-shaped optical 
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FIG. 3. The population imbalance in two coupled 

rings. We focused on the case A » Xp. For 

moderate zq, oscillations are obtained, with cj 
2g {vTTA2 + Xp(zoA - vT^)(2A2 - 1)/[2(1 + A^f^^]] 
corresponding to macroscopic quantum self trapping. The dynamics 
can be visualized with the help of the mechanical system provided 
by a rotator of length y^l - z{s)'^, driven by the external force 
A. The constant solution z{t) = const corresponds to vanishing 
pendulum length. For A = (inset), the dynamics is characterized 
by Rabi oscillation with uq = 2g{l + Ap\/1 - zo/2) < cj. Here 
Ap = 0.1 and A = 4 implying that uj « Aojq. 
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FIG. 4. The time of flight expansion for the two-coupled-rings-qubit. 
Left panel: vanishing inter-ring tunneling rate g/t = 0. Right panel: 
g/t = 0.9. In ihQ (kx^ky) plane the interference fringes with the 
ring symmetry are due to the momenta of the quantum degenerate 
gas; the inter-ring tunneling suppresses the interference fringes. In 
the (ky^kz) plane g, induces structured interference fringes. The Eq. 
(8) is calculated for the Bose-Hubbard ladder with 'fluxes' $« and 
$6, with [7 = 0, and at quantum degeneracy. In the plots: $a = 
80, $b = 70, T = OmkB and iV = 14 with filling fractions of 10 
bosons per site. 



lattice potentials. Persistent currents had been experimentally 
observed [37] in a narrow toroidal trap with a weak link. The 
effective action of the system studied in [38] can provide a two 
level system. In contrast with [37, 38], we emphasize how we 
make explicit use of the lattice in our construction, both to 
confine the particles in the rings and to drive the ring-ring in- 
teraction. The qubits are realized by breaking the Galilean 
invariance of the system either by adding an additional bar- 
rier along a single ring lattice Eqs.(2), or by tunnel coupling 
of two homogeneous rings, Eq.(5). The latter is proposed to 
be realized with the scheme in Fig. 2. At the degeneracy point 
basis, the effective Hamiltonians (2) and (5) lead to one and 
two qubit operations (see [39]). 

Our work provides a feasible route to the implementation of 
a functional flux qubit based on persistent atomic currents. 
The initialization of our qubit can be accomplished, for exam- 
ple, imparting rotation by exploiting light induced torque from 
Laguerre-Gauss (LG) beams carrying optical angular momen- 
tum. A two-photon Raman transition between internal atomic 
states can then be used to transfer coherently h orbital angu- 
lar momentum to the atoms. With this method, transfer effi- 
ciencies of 90% to the rotating state had been demonstrated 
[37, 40]. Owing to the coherent nature of the Raman process, 
superpositions of different angular momentum states can be 
prepared [41, 42]. Measurements of the decay dynamics of a 
rotating condensate in an optical ring trap showed remarkable 
long lifetimes of the quantized flow states on the order of tens 
of seconds even for high angular momentum (1=10)[41]. Con- 
densate fragmentation and collective excitations which would 
destroy the quantum state had not been observed. To allow 
controlled tunneling between neighboring lattice stacks the 
distance between the ring potentials needs to be adjustable in 
the optical wavelength regime (the schematics in Fig. 2 can be 
employed). Small distances allow high tunneling rates, a ne- 
cessity for fast gate operations. This makes it however less ef- 
ficient to read out and address individual stack sites. Increas- 
ing the lattice stack separation after the tunneling interaction 
has occurred well above the diffraction limit while keeping the 
atoms confined, optical detection and addressing of individual 
rings becomes possible. 

Read out of the angular momentum states can be accom- 
plished experimentally with interference of different flow 
states (i.e. corresponding to a fragmented superfluid) which 
maps the phase winding into a density modulation that can 
be measured using time-of-flight imaging [41]. In the lower 
panel of Fig.4 it is shown that different flow states lead to 
characteristic density patterns in the far field. 

We believe that our implementation combines the advan- 
tages of neutral cold atoms and solid state Josephson junc- 
tion based flux qubits for applications in quantum simulation 
and computation. This promises to exploit the typically low 
decoherence rates of the cold atom systems, overcoming the 
single site addressing, and harness the full power of macro- 
scopic quantum phenomena in topologically non trivial sys- 
tems. The characteristic fluctuations in the magnetic fields 
affecting Josephson junction based flux qubits are expected to 



be minimized employing neutral atoms as flux carriers. 
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